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Related AppiicatiQns 

[0001] This is a continuation application under 35 U.S.C. S 120 of WO 
2005/041449 AK filed as PCT/BE2004/00Q154 on October 27, 2004, which is hereby 
incorporated by reference. This application claims the benefit of the filing date of U.S. 
Provisional Patent Application No. 60/514,625, filed October 27, 2003. which is hereby 
incorporated by reference. 

Field of the invention 

{OOO44-[0002] The present invention is related to a method for charact e rising 
characterizing an optical fifere-fiber link by analysis of a Rayleigh backscatter signal, from 
which the spatial distribution of the beat length, the coupling length and the PMD 
(Polarization Mode Dispersion) along the link can be determined. 

State of the art 

{0003}-[0003] The use of dispersion-shifted and dispersion compensating fibres 
fibers has minimised the effect of chromatic dispersion on the bandwidth of an optical 
link. Polarisation Polarization mode dispersion (PMD) has therefore become the most 
serious limiting factor in high-speed optical communication systems. There exist several 
measurement techniques of PMD, but they only allow for measuring the global value and 
do not give information about the distribution of PMD along the fibfe-fiber length. 
Measuring the spatial distribution of PMD is, however, important. Such a measurement 
would indeed allow to locate the bad trunks in an optical link, causing high PMD value. 
This can be essential in the frame of a network maintenance and upgrade. 

{000^[0004) The PMD of a fifefe -fiber depends on two parameters: the beat 
length (Lb) and the coupling length (Lc). The beat length depends on the birefringence of 
the fibre -fiber and the coupling length is related to the phenomenon of polarisation 
polarization mode coupling. Measuring the spatial distribution of PMD therefore consists 
ef involves measuring the distribution of both the beat and coupling lengths. 
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f0004f[00051 Most of the techniques proposed for measurement of the distributed 
PMD are based on polari s ation polarizat ion -opt icai time domain reflectometry (POTDR). 
The concept of POTDR was introduced by A. Rogers in the 80's (see A. Rogers, 
'Polarization-optical time domain reflectometry: a technique for the measurement of field 
distributions'. Applied Optics, vol. 20, pp. 1060-1074, 1981). It basically consists of 
involves measuring the polarisation polarization properties of the Rayleigh backscatter 
signal when an optical pulse propagates down the fibre -fiber . Some POTDR-based 
techniques have since been developed for PMD characterisation characterization . Most of 
them allow a measurement of the beat length distribution along an optical link and, 
therefore, do not allow the complete determination of the PMD. However, Sunnerud 
presented a POTDR-based set-up enabling the measurement of the accumulation of PMD 
along the fibre — fiber (see H. Sunnerud et ai, Tolarization-Mode Dispersion 
measurements along installed optical fibers using gated backscattered light and a 
polarimeter\ 7. of Lightwave TechnoL, vol 18, pp. 897-904, 2000). This method requires 
POTDR measurements at several wavelengths. A technique based on the measurement of 
the degree of polarisation polarization by using a POTDR set-up has also been presented 
in 'Distributed detection of high-PMD sections on installed fibers using a polarization- 
OTDR\ M. Leblanc, Proc, OFMC'OL pp. J 55-1 62, 2001, and is related to the patent 
application US 2003/0174312. This technique does not allow quantification of the PMD, 
but enables to obtain distributed information about the level of PMD. Moreover this 
technique requires the complete measurement of the state of polarisation polarization of 
the backscattered signal. More recently, Galtarossa described another technique for the 
PMD measurement for which the coupling length is determined from the correlation of 
the birefringence vector (see 'Measurement of birefringence correlation length in long, 
single-mode fibers', Galtarossa et aL, Opt. Lett,, vol. 26, pp. 962-964, 2001), which also 
requires the measurement of the polarisation polarization states of the backscattered 
signal. Moreover, this technique requires a quite small duration of the pulses launched in 
the optical fibre -fiber : 5-10 ns. 

[0005] [0006] In 'Polarization mode dispersion mapping in optical fibers with a 
Polarization-OTDR M. Wuilpart, G. Ravet, P. Megret and M Blondel, Photon. Technol. 
Lett., vol. 14, pp. 1716-1718, 2002, the authors describe an analysis of the POTDR signal 
which allows a mapping of the PMD along an optical fibre -fiber link by quantifying the 
PMD on each fibre-fiber of the link. The main advantage of that technique is that it does 
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not require the complete measurement of the backscattered polariGation polarization state 
evolution and only uses a linear polariser polarizer at the fibfe-fiber input, which is quite 
simple to implement. The determination of the beat length, the coupling length and finally 
the PMD is based on the analysis of the statistical properties of the extrema present in a 
POTDR trace. 

{0006]-[0007] In said paper by Wuilpart et al. the spatial distribution of the linear 
birefringence 5 is characteris e d characterized by a Rayleigh statistical distribution along 
the fibfe-fiber length and it is assumed that the rate of change of the birefringence angle q 
is described by a Gaussian distribution with zero mean and standard deviation a 
depending on the polarisation polarization mode coupling strength. The circular 
birefringence is neglected. The corresponding measurement set-up is shown in Fig. 1, The 
OTDR JTO pulses (10 ns) externally modulate a 1553 nm distributed feedback ( DFB) 
laser 120 via a pulse generator 130 and an electro-optic modulator (EOM) 140 . After 
amplification, the pulses are launched into the fibr e fiber through an acousto-optic 
modulator (AOM) 150 , which suppresses the amplified spontaneous emission noise of the 
erbium doped fifefe -fiber amplifier between two successive pulses. A linear polaris e r 
polarizer 160 is placed at the fibfe -fiber input and a polarisation polarization controller 
(CP) 170 is used to obtain the maximum power after the polariser polarizer . 

{000^[0008] Fig. 2 shows three simulated POTDR traces corresponding to the 
POTDR arrangement and obtained for different values of the beat length and an identical 
a value. One can easily observe that the number of maxima strongly depends on Lb. 

|#008j-[00091 The calculation of the beat length is based on the measurement of 
the number of maxima of the POTDR trace. Fig. 3, which has been obtained by 
simulation, presents the number n of maxima per unit of length in function of the beat 
length for different values of a. This figure shows that n varies not only with Lb but also 
with the strength of the polarisafion polarization mode coupling (a). In order to estimate 
Lb from n, the curve corresponding to a = 10 (the dashed curve in Fig. 3) is used. This 
choice indeed minimises the maximum possible error if one assumes that a is smaller 
than 15 degrees/m. Consequently, this method allows only to give an approximate value 
of the beat length Lb. A relationship between the coupling length Lc and a can be derived 
and one can obtain: 
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{00&9f [0010] A statistical approach is then used for determining the coupling 
length from the POTDR trace. Fig. 4 shows different POTDR traces that were obtained 
by simulation for a beat length of 25 m and for different values of a: 2, 4 and 8 
degrees/m, respectively. One can observe that the lower envelope of the curves varies in a 
smoother way when a decreases and therefore when Lc increases. Defining APi as the 
absolute value of the difference in ordinate between two successive minima of the 
POTDR trace P one can write : 

AP,=|/'(z,)-P(z,„)| 

where z\ and Zj+i correspond to the distances of the i and the i+l minimum, 
respectively. 

{0040]-[0011] A new parameter ^ can then be calculated as 

4 =< AP, >, 

where <.>i denotes the mean value for i varying from 1 to the number of minima of the 
POTDR trace. ^ is a measure of the speed variation of the lower envelope of the 
backscattered signal along the fibfe -fiber length. Fig, 5 (obtained by simulation) shows 
the evolution of J| according to a for different values of the beat length Lb. The values are 
here expressed in dB because £, has been calculated directly from the POTDR trace. 

{M14j-[0012] The procedure for the determination of Lc is the following. From 
the POTDR trace, the number of maxima and each APi are measured. An approximate 
value of the beat length is then deduced by the method previously described and ^ is 
calculated by averaging APj. Finally a can be deduced from )| by using the curve of Fig. 5 
corresponding to the correct beat length and then Lc is calculated by means of the relation 
Lc=0.879/a^ stated above. Because Lb is an approximate value, the calculated coupling 
length Lc will also be an approximation. 

f00i2}-[0013] From the measurement of the beat length and the coupling length, it 
is easy to determine an approximate value of the PMD value by 



PMD^ = i 
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where c is the light velocity in a vacuum, \ the measurement wavelength and L the fibre 
fiber length. It should be emphasised that the final PMD value is approximate: this is due 
to the approximation done on the beat length measurement. This process can be repeated 
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for each contribution in the POTDR signal corresponding to the different -fibfes -fibers of 
the optical link: it is therefore possible to map the PMD along this link. 

{M13]-(0014] The method as described above has some important drawbacks and 
limitations. It only yields approximate values of Lb, Lc and the PMD. Further, the 
presented POTDR traces were assumed ideal. It supposes ideal characteristics of the 
various components of the experimental set-up. in practice, the set-up is charact e ris e d 
characterized by some imperfections, which cause a distortion of the ideal POTDR signal. 
Therefore the method as such cannot be directly applied. The POTDR signal is also 
affected by noise: detector noise and residual coherence noise, which adds a series of 
minima and maxima to the POTDR trace. The final value of the measured PMD can 
therefore be erroneous. 

Aims of the invention -Summary of Certain Embodiments 

fOW4]-[0015] Th e pr e s e nt inv e ntion aims to Certain embodiments provide a 
method for analysing a backscatter signal for determining in an accurate way the spatial 
distribution of the beat length, the coupling length and the polarisation polarization mode 
dispersion. 

Summary of th e inv e ntion 

{OOiSJ-[0016] Th e inv e ntion relat e s Certain embodiments relate to a method fer-of 
characterising characterizing an optical fibre- fiber link by its beat length, coupling length 
and polarisation polarization mode dispersion distribution^ — It comprises the steps of 
comprising 

- sending a pulsed signal along the optical fibre fiber link and measuring the 
backscattered signal (preferably a POTDR signal), after passing through a polari se r 

polarizer, 

- deriving the length of the optical fibre-fiber, the average power difference between 
two successive minima of the backscattered signal and the number of maxima per unit 
length, 

- in an iterative way determining a beat length interval and an interval for the 
polarisation polarization mode coupling parameter, until the length of the intervals is 
below a predetermined value, yielding a value for the beat length and the coupling 
length, and 
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- calculating the polarisation polarization mode dispersion. 

tOW4}-[0017] Preferably the POTDR signal is an ideal POTDR signal. 

{0©1^[0018] Advantageously the POTDR signal is the convolution of an ideal 
POTDR signal and a signal depending on the pulse shape. In an alternative embodiment^ 
the POTDR signal further is convoluted with a signal taking into account the time jitter 
introduced by the measurement set-up. 

{OW8]-[0019] In a specific embodiment a smoothing algorithm is applied to the 
POTDR signal. 

fOW-9j-[0020] In another obj e ct the inv e ntion Another embodiment relates to a 
method for characteri s ing characterizing an optical link consi s ting of comprising a 
concatenation of several fibr e s fibers, wherein the method as previously described is 
applied to each fibre-fiber. 

{0020]-[0021] The inv e ntion also Another embodiment relates to the use of a 
method as described above to locate the position of polarisation polarization mode 
dispersion sources within an optical fibre -fiber link. 

{002ij-[00221 In an advantageous embodiment^ the method as described is used in 
telecommunication networks. Alternatively, it can be applied in fibre fiber sensing 
applications. 

Short description of the drawings 

{00221-[0023] Fig. 1 repr e sents the is a block diagram of an example set-up for 
measuring a POTDR trace. 

{002^(0024] Fig. 2 repres e nts is a diagram of three simulated POTDR traces for 
a given a and for different Lb values. 

{0024]-[0025] Fig. 3 r e pr e s e nts is a diagram of the number of maxima per unit 
length versus Lb for different values of a. 

{002§}-[0026] Fig. 4 represents is a diagram of simulated traces for a given beat 
length and different values of a. 

{0§2^[0027] Fig. 5 repr e s e nt s is a diagram of the evolution of as a function of 
a for different values of Lb. 

{002^(0028] Fig. 6 r e pr e s e nts is a diagram of the determination of the beat 
length interval. 



6 



{002S]-[0029] Fig. 7 represents is a diagram of the determination of the interval 

for a. 

{002^[0030] Fig. 8 repr e s e nts is a diagram of the distorsion distortion of the 
POTDR signal due to the imperfections of the set-up (pulse shape -i- time jitter). 

{OOSO]-[0031] Fig. 9 represents is a diagram of the evolution of n with the beat 
length when the imperfections of the set-up are taken into account. 

{003if|0032] Fig. 10 represent s is a diagram of the evolution of 4 with a when 
the imperfections of the set-up are taken into account. 

{W33j-[0033] Fig. 1 1 represents is a diagram of the simulated and measured 
probability density function of A. 

[0033] [0034] Fig. 12 r e presents is a diagram representing the window algorithm 
that can be used for calculated the local PMD. 

Detailed description of the- kiv e ntion certain embodiments 

{^034]-[0035] It is first shown how the polarisation polarization mode dispersion 
from a POTDR signal can be determined in an accurate way (Lb is not approximately 
determined anymore) by analysing the statistical properties of its extrema. Suppose first 
that a falls within the range [0,15] degrees/m, which means that the coupling length is 
assumed to be larger than 12.8 m. This a interval has been chosen such that it 
corresponds to the values measured on real fibres -fibers . 

{00351-[0036] From a POTDR signal, the fibre-fiber length L, ^ and the number of 
maxima per unit of length, nm, are easily measurable. After measuring nm, it is possible to 



determine a first interval 



that: 
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and are determined by the intersections of the curves n(LB,a) corresponding 

to a = 0 and a = 15 degrees/m with the straight line n = nm as shown on Fig. 6. 

{003^(0037] When the first beat length interval is determined, the curves 

^(o,Lb) can be used to define an interval P'/ow^^ip. for the a value of the fibre -fiber 
under test so that : 
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o'iow and <^lp are determined by the intersections of the curves ^(a^Le) for - 
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and ^is " with the straight line ^ = ^m, where is the ^ value of the fibre -fiber 

up 

under test measured from the POTDR signal. This principle is illustrated in Fig. 7. 

f003^[0038] After the determination of (^hw and crip, the curves n(LB,a) for 

^ ~ ^Lw and = can be calculated. A new beat length interval is then deduced as 

previously explained by Fig. 6, but using the curves obtained for o" = ct/^^ and cr = 
instead of a = 0 and a = 15 degrees/m. 



f0038}-[0039] This new beat length interval denoted 
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and gives a smaller range for the beat length of the fibre -fiber under test. 



Fig. 7 is then used again to calculate a new a interval P/wCr;;^] by using the curves 



i^(a,LB) obtained for 



^B - ^ o and ^B - ^ o 
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This new interval gives a smaller range for the a value of the fibre -fiber under test. This 
process is repeated until small enough intervals are obtained for the beat length Lb and for 
a. 

fOO344-[00401 When the beat length Lb and a are determined, the PMD is finally 
calculated by means of the relationship: 
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f0040}-[0041] If the optical link to be characteri se d c haracterized 

consists of comprises a concatenation of several fibres-fibers, the method is applied to 
each fibre -fiber of the link and their beat and coupling lengths are therefore determined 
independently of each other. In this way, the PMD mapping is possible along the link. 

[OO44-|-|0042J The method for analysing the backscatter signal as 

described offers several advantages : 
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- it requires only the access to one end of the optical fibfe- fiber link to be charact e ris e d 
characterized : 

it is based on a POTDR arrangement, but does not require the complete measurement 
of the complete state of polarisation polarization of the backscattered signal. Only a 
simple polariser polarizer is required; 
it involves a POTDR measurement at only one wavelength; 

- it gives accurate (not approximate) values for the beat length, the coupling length and 
the PMD. 

f0043}-[0043] In practice, the final measured signal is also affected by noise. This 
noise is due to the detector noise and to the presence of a residual coherence noise. As the 
determination of the PMD is based on the detection of extrema, it is obvious that this 
noise will affect the results. 

{004^(0044] As already mentioned before, the measured signal is also affected 
by the imperfections of the experimental set-up like the finite pulse shape and the time 
jitter of the electronic equipment. 

100451 These phenomena (noise and imperfections of the measurement set-up) 
lead to a distortion of the POTDR signal. Therefore it cannot be used directly to deduce 
the PMD map. 

100461 The width of the optical pulses launched into the fibre -fiber under test 
is finite and can be non-rectangular. Consequently, the backscattered power measured at 
the POTDR corresponds to an average of the ideal signal over a small spatial window 
depending on the pulse shape. The resulting signal measured at the POTDR is in fact the 
convolution product between the ideal signal and the shape of the pulse: 

where Ppuise(z) is the power measured by the POTDR according to a certain pulse shape. 
P(z) is the ideal POTDR signal and Pshape(z) is the power distribution along the pulse. 
Hence the Fourier transform of Ppuise(z) becomes: 

P pulse 

where Ppuise(k), P(k) and Pshape(k) are the Fourier transforms of PpuiseC^), P(z) and Pshape(z), 
respectively. 

{0044J-(0047] Moreover the electronic equipment is not perfect: it is subject to a 
time jitter, iter e.g., the time delay between two successive pulses is not constant. This 
phenomenon strongly affects the averaging process of the POTDR. If we suppose that the 
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number of pulses involved in the averaging process of the OTDR is large enough and that 
the time jitter is a random variable characterised characterized by a Gaussian distribution 
of zero mean and standard deviation Oy, the signal detected at the POTDR can be written 
as: 

where Pm(z) is the power measured by the POTDR taking into account the pulse shape 
and the time jitter. Hence the Fourier transform of Pm(z) becomes: 

The pulse shape and the time jitter have for effect to distort the measured backscattered 
signal. An example is shown on Fig. 8, which has been obtained for parameters of the 
experimental set-up. The signals (a) and (b) correspond to the ideal and affected signals, 
respectively. One can clearly observe that the signal is strongly distorted. 

{004S}-[0048] Because of these imperfections, the method described above cannot 

be used directly: the curves n(LB,a) and ^(OjLb) are indeed based on ideal POTDR 
signals obtained by simulations. Therefore the effects of the pulse shape and the time 
jitter should be included in the calculation of n(LB,a) and ^(a,LB). In the case of the 
parameters of the experimental set-up and a pulse duration of 10 ns, these new curves are 
shown on Fig. 9 and 10. In order to deduce these curves, ideal POTDR traces have been 
simulated and the effects of the imperfections have been taken into account by using the 
above expressions for Ppuhe(z) and Pm(z)- Consequently, when the imperfections of the 
experimental set-up are taken into account, the PMD map of an optical link can still be 
determined. The same reasoning as the one explained previously has to be applied on the 
POTDR signal by using the new set of curves n(LB,a) and E,(a,LB) described in figr -Figs. 
9 and 10. As the pulse shape is now taken into account, it is possible to launch a larger 
pulse in the fibre -fiber (50 ns, for example), which will increase the system dynamics. 
The maximum measurable length will therefore be greater. 

f0046}-[0049] As already indicated, the final measured signal is also affected by 
noise. This noise is due to the detector noise and the presence of a residual coherence 
noise. Its main effect is to add a series of minima and maxima to the POTDR trace. The 
values of n and 4 cannot thus be correctly measured and the final calculated PMD will be 
erroneous. In order to solve this problem, a smoothing algorithm is applied to the 
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backscattered POTDR signal Pni(z). This algorithm consists of comprises generating the 
signal Ps(z) such that: 



Ps(2) is therefore the mean value of Pm(z) on a distance interval Ls around z. The number 
of maxima and the distances Zk corresponding to the locations of the minima of the 
POTDR trace are then calculated from the smoothed signal in order to ignore the 
influence of the small peaks generated by the noise. Afterwards, when the positions Zk are 
known, ^ is calculated from Pm(z). When noise is present, the number of extrema of Ps(z) 
depends of the smoothing parameter Lj. So, the main problem is to choose the correct 
smoothing parameter to be applied to the POTDR signal. 

{0047J-[0050] The location of the minima of the smoothed signal can be used to 
estimate the probability density function of A, denoting the difference in z between two 
successive minima. This estimation can then be compared to the density probability 
function of A of a POTDR trace obtained by means of simulations taking into account the 
parameters Lb and a calculated from the smoothed signal. If the probability density 
functions of A obtained from Ps(z) and from the simulation do not match, it means that 
the noise contribution is significant. The smoothing parameter Ls will therefore be chosen 
such that the corresponding probability density function obtained from measurements can 
be correctly fitted by the probability density function derived from simulations: Ls is 
therefore chosen such that the noise contribution is reduced. Fig. 1 1 illustrates this 
concept. Figures ii(a), ii(b) and ii(c) show the probability density functions obtained 
by measurement (circles) and by simulation (crosses) for a given optical fibre -fiber and 
for an Ls value of 3, 4 and 5 m, respectively. One can clearly observe in this case Ls = 3 m 
gives the best fitting. 

f0048f [0051] Th e inv e ntion also provid e s Certain embodiments provide a 
technique which will allow allows to locale more accurately the position of the PMD 
source within the fifere-fiber. The basic idea consists in involves applying the described 
method on a distance window (1 km, for example), which slides along the fibre-fiber 
length. By comparing the obtained PMD's on several successive windows as shown on 
Fig. 12, it is possible to locate more accurately high PMD sources within the optical fibre 
fiber. 
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t0049f [0052] POTDR can also be used in fibre -Fiber sensing. The polarisation 
polarization properties indeed depend on several external effects like strain, temperature, 
electric and magnetic fields. Strain, twist, temperature, electric and magnetic fields affect 
the birefringence within the fibre-fiber and consequently, the mean beat length. The 
method of th e inv e ntion Certain embodiments, which measures measure the spatial 
distribution of the beat length, therefore allows allow to measure the spatial distribution 
of these external effects after calibrating the sensing fibre-fiber. 
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